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Abstract 

The fourfold degeneracy of the boron acceptor ground state in silicon, which is easily lifted 
by any symmetry breaking perturbation, allows for a strong inhomogeneous broadening of the 
boron-related electron paramagnetic resonance (EPR) lines, e.g. by a random distribution of local 
strains. However, since EPR of boron acceptors in externally unstrained silicon was reported for 
the first time, neither the line shape nor the magnitude of the residual broadening observed in 
samples with high crystalline purity were compatible with the low concentrations of carbon and 
oxygen point defects, being the predominant source of random local strain. Adapting a theoretical 
model which has been applied to understand the acceptor ground state splitting in the absence of 
a magnetic field as an effect due to the presence of different silicon isotopes, we show that local 
fluctuations of the valence band edge due to different isotopic configurations in the vicinity of the 
boron acceptors can quantitatively account for all inhomogeneous broadening effects in high purity 
Si with a natural isotope composition. Our calculations show that such an isotopic perturbation 
also leads to a shift in the g-value of different boron-related resonances, which we could verify in 
our experiments. Further, our results provide an independent test and verification of the valence 
band offsets between the different Si isotopes determined in previous works. 
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I. INTRODUCTION 



Electron paramagnetic resonance (EPR) and related techniques like electron nuclear dou- 
ble resonance have contributed extensively to the understanding of substitutional shallow 
donors in the different allotropes of silicon throughout the last 50 years.- - - On the contrary, 
EPR had been ineffective for the study of shallow acceptors for a long time. The reason 
behind this asymmetry lies in the different structure of the respective dopant ground states. 
While the electronic ground state of substitutional shallow donors in silicon is s-like and 
only twofold spin degenerate, shallow acceptors states have p-character and show a four- 
fold degeneracy- This latter degeneracy can partially be lifted by any symmetry breaking 
perturbation. Such perturbations, e.g. electric fields or strain, can strongly alter or even 
dominate the level scheme of acceptor Zeeman energies for external magnetic fields that 
are typically used for EPR measurements. If the perturbation is not homogeneous across 
the sample, this can easily lead to an extreme inhomogeneous broadening of the EPR res- 
onances. Therefore, the observation of a boron-related EPR signal had initially only been 
possible under application of a strong and homogeneous external stress.- 

It was only in 1978, when Si became available with sufficient crystalline quality, that 
EPR of B acceptors in externally unstrained Si (Si:B) was reported for the first time by 
Neubrand.— ^ In particular, a correlation between the linewidths of the different boron- 
related EPR resonances and the concentrations of C and O point defects was established. 
However, a number of fundamental questions have remained unsolved: (i) Although a ran- 
dom strain distribution that is induced by point defects should lead to a purely Lorentzian 
broadening of the EPR resonances, the experimentally obtained lines could only be fitted 
with Voigt profiles, taking into account a significant Gaussian contribution that was found to 
be independent of the concentrations of C and O. The origin of this additional broadening, 
which dominates the overall linewidths in samples with small point defect concentrations 
and which shows a large angular dependence, has essentially remained unclear, (ii) For the 
samples of highest crystalline purity, the Lorentzian contribution to the Voigt profiles did 
not fall below a threshold value of 10 mT which, in the model developed by Neubrand, 
would correspond to an unreasonably high concentration of point defects. In Ref. [ill, Si 
interstitials were suggested as a possible explanation for this finding, (iii) For two of the 
six B-related resonances, a distinct substructure was observed. It was proposed that this 
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structure originates from a dynamic effect, however, the specific mechanism remained un- 
known.—*^ 

Karaiskaij and coworkers proposed in Ref. 
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that the random distribution of the different 



stable Si isotopes 28 Si, 29 Si, and 30 Si in Si crystals with a natural isotope composition 
plays an important role for the inhomogeneous broadening of B-related EPR resonances. In 
that work, it was shown that the valence band offsets between isotopically pure Si crystals 
consisting of different isotopes, which lead to local fluctuations of the valence band edge in 
the vicinity of the different B acceptors in nat Si, are responsible for the residual ground state 
splitting of shallow acceptors in this material, which is e.g. observed in photoluminescence 
spectra of acceptor-bound excitons.— ^ However, it has remained unclear which of the open 
questions listed above can be attributed to local valence band fluctuations and whether a 
quantitative understanding of the observed effects is possible. In the present paper, we 
extend the theoretical model established in Ref. ljj to explain the acceptor ground state 
splitting in the absence of an external magnetic field (B=0) to non-zero magnetic fields and 
investigate the influence of isotope-induced perturbations on the line shape of B-related EPR 
resonances in Si with the natural and with isotopically engineered isotope compositions. We 
show that the inhomogeneous broadening effects described in (i) and (ii) can quantitatively 
be explained by isotope-induced random local fluctuations of the valence band edge in the 
vicinity of the B acceptors. We verify this finding experimentally via a direct comparison 
of EPR spectra obtained from B-doped nat Si and isotopically purified B-doped 28 Si. The 
isotope-induced broadening is discussed for different orientations of the magnetic field and 
different isotope compositions. Additionally, the asymmetric line shape of the B-related EPR 
resonances observed in nat Si, which was not addressed in previous work can quantitatively 
be explained using our model. We show that the random distribution of the different Si 
isotopes in nat Si also shifts the effective g-values of different boron-related EPR resonances 
with respect to pure 28 Si. This effect is quantitatively verified by our experimental data. 
The crucial parameters for the agreement between this model and the experimental data 
are the valence band offsets between pure 28 Si, 29 Si, and 30 Si, which were extracted from 
calculations of the temperature dependence of electronic band states.— 1 ^ Using our model 
to recalculate the residual acceptor ground state splitting for B=0, we show that the 30 % 
discrepancy between theory and the experimental value determined by phonon absorption 
spectroscopy in Ref. [l4| is not due to an uncertainty in the values assumed for the valence 
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band offsets, but rather due to an inaccuracy of the acceptor wave function used in that 
publication. In our calculations, the best agreement with the experimental data is obtained 
when the assumed valence band offset between 28 Si and 29 Si is reduced by 8 % from the 



II. THEORETICAL MODEL 

Boron acceptors in Si are conventionally described in the hydrogenic impurity model, 
where a hole is bound by the electrostatic potential of the singly charged B-ion which 
is screened by the dielectric constant of Si.- The ground state envelope function of the 
hole is the Is state of a hydrogen-like atom. As shown in Fig. [T]^a) for B in pure 28 Si, 
this ground state is fourfold degenerate at B=0 due to the spin and Bloch basis angular 
momentum.- This degeneracy is partially lifted in nat Si due to the action of isotope- induced 
local fluctuations of the valence band edge in the spatial vicinity of the different B acceptor 
nuclei which constitute a symmetry breaking perturbation for the acceptor wave function. 
These fluctuations result from the dependence of the zero-point vibrational energy on the 
isotope mass in conjunction with the renormalization of the electronic states by electron- 
phonon interaction.— In the ensemble average, the resulting two Kramers degenerate levels 
are shifted by an energy that equals the shift of the valence band edge from 28 Si to nat Si. 
However, for different B acceptors, the shift of the different energy levels will undergo slight 
variations depending on the specific distribution of Si isotopes in the vicinity of the acceptor 
nucleus. 

When a strong external magnetic field is applied, the situation is different, as shown 
in Fig. []Jb). Already the Zeeman interaction between the acceptor-bound holes and the 
external magnetic field leads to a full lifting of the degeneracy of the heavy and light hole 
states and the energy eigenstates can approximately be classified in the total angular mo- 
mentum basis \j=3/2,m e {±1/2, ±3/2}), where j and m are the angular momentum and 
the magnetic quantum number, respectively. For strong enough external magnetic fields, the 
isotope-induced effects can then be treated as a small perturbation. As schematically shown 
in Fig. Hfb), on average this perturbation should lead to an upwards shift of the four Zeeman 
levels in energy similar to the situation shown in Fig. Ufa) for B=0. The assumption that 
this picture holds for X-band EPR conditions is motivated by the experimental evidence 



value determined in Refs. 14. 
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FIG. 1: (Color online only) (a,b) Energy level schemes of the acceptor ground state, (a) In the 
absence of an external magnetic field B, the fourfold degenerate ground state in 28 Si is split into 
two twofold degenerate states in Si both shifted to higher energies, (b) When a magnetic field 
of e.g. B = 0.6 T is applied, the degeneracy of the ground state is already fully lifted by the 
Zeeman coupling in 28 Si. The isotopic effect leads to an additional upwards shift of the energy 
levels which, for a given isotopic configuration, slightly varies for the different Zeeman levels, (c) 
Two-dimensional cross sections at z=0 (acceptor situated at the origin) of the envelope functions 
of the acceptor ground state probability for the four Zeeman states shown on a logarithmic scale. 



that the measured electronic g- values of B in nat Si are predominantly determined by the 
Zeeman interaction in samples with a low concentration of point defects.— The isotopic dis- 
order therefore only leads to a comparatively small perturbation that mostly manifests itself 
in a change of the linewidths and lineshapes of the different boron-related EPR resonances. 

To assess the nature and the magnitude of these isotope-induced changes, we implement 
our theoretical model in the following way: First, we calculate the four Zeeman levels of 
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the acceptor ground state using a k • p envelope function model taking into account all six 
(2x3) spin- resolved valence bands including the split-off band. The Hamiltonian of the 
acceptor-bound hole can schematically be written in the form 

H = + ^s 6x6 -B + V(r), (1) 

which was solved on a real-space grid using the NextNano++ code.— The first term on the 
right hand side represents the six-band effective mass Hamiltonian in a discrete real-space 
basis including the magnetic field B in a non-perturbative and gauge-invariant manner, with 
B only appearing in phase factors. The second term, where hb is the Bohr magneton, g is 
the free electron g-value and the 6x6 spin matrices Sk=l 3x3 ®ak(ke{x, y, z}) are determined 
by the Pauli matrices &k, couples the spin to the magnetic field. The impurity nucleus is 
represented by a negative charge at the center of the simulation domain screened by the bulk 
silicon dielectric constant of £si=H-7. The potential energy of the acceptor hole is given by 

e 2 

V{r) = + W(r), (2) 

where we have included the so-called central cell correction W(r). It phenomenologically 
corrects for deviations of V(r) from a purely hydrogenic potential close to the impurity 
nucleus due to changes in the dielectric properties in this spatial region.^ To account for 
W(r), we adopt the parametrization 



e 2 



W(r) = — 



Ae { ~ ar) + (1 - A)e 



c 



(-7T-) 



(3) 



used in Ref. [18|, where the parameters a=0.755/a B , /3=0.35/aB, 7=2.45/a B , and v4=1.14 
are phenomenological fitting parameters and a-Q is the Bohr radius of the effective mass 
acceptor.— In order to avoid a singularity of V(r) at r=0, the Coulomb potential is replaced 
by QS rj o for the central grid node following Ref. [l8|. The only free parameter in our model, 
Q, has been chosen to reproduce the correct mixing of the different valence band states in 
which is responsible for the inhomogeneous energy splitting between the Zeeman levels 
in the absence of any further perturbations. A direct experimental measure for this mixing is 
the difference between the transition energies of the transitions l^i)*-?-! 1 !^) and l^)^!^)- 
This energy difference can directly be calculated from the difference in the resonance fields 
of the broad and the narrow Am=l EPR lines observed in 28 Si. We have fitted Q to 
reproduce the experimentally observed splitting of 49 mT for B \ \ [001] (cf. Fig.[3^b), Sec. lIV)) . 



For all further calculations, Q has been fixed to this value. In order to account for the 
long range character of V(r) as well as the strong confinement at the acceptor nucleus, we 
use a simulation domain size of 50x50x50 nm 3 and an inhomogeneous grid with a strong 
concentration of nodes close to the acceptor position, respectively. The wave functions 
were forced to zero at the boundaries of the simulation domain. With this model, we 
have calculated the four ground state Zeeman levels labeled by 1^2), 1^3), an d l^) 

(corresponding to m=— 3/2, m=— 1/2, m=+l/2, and m=+3/2, respectively), as shown in 
Fig. Hfb), for external magnetic fields with different field strengths and orientations relative 
to the crystallographic axes. In Fig. (He), we exemplarily show typical wave functions 
obtained by these calculations for B=0.6 T with 5=5 2 ||[001] as two-dimensional cross 
sections through the different Zeeman states at z=0 and the B nucleus located at the origin. 
The band warping of the different valence bands is reflected in the structure of the envelope 
functions. While the two heavy hole- and light hole-like states are very similar to each 
other, between them there are distinct differences. From this it already becomes clear that 
the zeroth order effect of any perturbation will be an energy shift of the inner two light 
hole-like states ^2) and ^3) with respect to the outer two heavy hole-like states j^i) and 
^4) as schematically indicated in Fig. U^b). 

Based on the calculated wave functions, the isotope-induced changes of the band edges 
are considered perturbatively. For this, we map the wave functions on a discrete atomistic 
silicon lattice. The values of the acceptor envelope functions ^(n) at the position of the 
n-th silicon site are evaluated by linear interpolation. For a silicon lattice that comprises 
42x42x42 unit cells, corresponding to the size of the simulation domain, we first choose 
random configurations of 28 Si, 29 Si, and 30 Si isotopes, which on average corresponds to the 
isotope composition intended for the different calculations, e.g. 92.23 % 28 Si, 4.67 % 29 Si, 
and 3.1 % 30 Si for natural Si. Following Ref. we introduce the isotopic perturbation 
potential 

for 28 Si 

VM = I AE 29 for 29 Si • (4) 
AE 30 for 30 Si 

The valence band offsets A£' 29 =0.74 meV and AE 30 =1A6 meV are taken relative to the 
28 Si valence band edge and were deduced from calculations of the temperature dependence 
of the electronic band states in Si.— >^ Next, we project the diagonal perturbation potential 
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into the subspace spanned by the four Zeeman states leading to a 4x4 perturbation 
Hamiltonian H?* 4 with the matrix elements 



Finally, we diagonalize the total Hamiltonian H to t = + H 4 S * 4 including the diagonal 

Zeeman term Hj^ = E^ij, with Ei being the energy of the z-th Zeeman level determined 
from the k • p model described above. The isotopic contribution H 4 ,* 4 leads to a mixing of 
the four Zeeman levels and to an upward shift of the states, as depicted in Fig. [Tfb). The 
relative alignment of the four resulting states |1), |2), |3), and |4) (cf. Fig. [TJb)) strongly 
depends on the specific isotopic configuration in the vicinity of the acceptor nucleus. We 
have therefore performed the calculation of the isotope shifts as described above for 200,000 
different random isotope configurations, taking into account different magnetic fields, isotope 
compositions, and isotopic perturbation potentials. The resulting statistical distributions 
of EPR transition energies allow the quantitative comparison with the experimental data 
presented below. 

III. EXPERIMENTAL DETAILS 

The nat Si:B sample studied has a B doping concentration of 10 14 cm -3 , was grown using 
the float-zone technique, has a rectangular shape and a size of approximately 3x3x9 mm 3 . 
The isotopically purified sample, referred to as 28 Si:B in the text, comes from the neck region 
of a float-zone crystal and has an enrichment of 99.98 % 28 Si. It is doped with boron to a 
concentration of 3(±l)xl0 14 cm -3 , has a cylindrical shape with a length of 10 mm in [100] 
direction and a diameter of ~3 mm. 

EPR measurements were performed with a Bruker Elexsys E 500 spectrometer in con- 
junction with a super-high-Q resonator (ER-4122SHQE) operated at an X-band microwave 
frequency of z/~9.4 GHz. During the measurement, the samples were cooled to temperatures 
of typically 3 K using an Oxford ESR 900 helium-flow cryostat. 
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FIG. 2: Statistical distribution of the transition energies between the different Zeeman levels for 
nat Si. The transitions are grouped according to the change in the magnetic quantum number Am. 
The dashed lines show the transition energies for pure 28 Si, as obtained from our k • p model 
without isotopic perturbation. The calculations were performed for an external magnetic field 
i?=0.6 T oriented along the [001] axis. 

IV. RESULTS AND DISCUSSION 



In Fig. |2]the statistical distributions that result from the calculations described above are 
shown for the six different EPR transitions in nat Si for B\\ [001]. We find that the distribution 
of the inner Am=l transition (|2)<H-|3)) and the Am=3 transition (|l)o|4)) are relatively 
sharp, while the outer Am=l (|1)-B-|2) and l 3 )^H 4 )) and the Am=2 (|1)<B-|3) and |2)<h-|4)) 
transitions show a full width at half maximum (FWHM) of approximately 10 % and 5 % of 
their mean value, respectively. This can be understood from the fact that the two former 
transitions occur between two states that are dominantly light holes or heavy holes, respec- 
tively, and experience a similar shift due to the isotopic perturbation. These transitions are 
therefore called intrasubband transitions. In contrast, the latter intersubband transitions 
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comprise both light hole-like and heavy hole-like states that are differently perturbed by the 
random isotopic distribution. The asymmetric broadening of the |1)-H-|4) and the |2)<B-|3) 
transitions results from the mixing of the four unperturbed eigenstates |^) that is obtained 
by the diagonalization of the total Hamiltonian H tot . As can be seen, the statistical dis- 
tribution of the two outer Am=l transitions and the Am=2 transitions are also slightly 
asymmetric and do not exactly overlap with each other. For low temperatures, when spin 
polarization becomes significant and the equilibrium populations of the four Zeeman states 
strongly deviate from each other, this result should be reflected in an increasing asymmetry 
of the EPR lineshape. The dashed lines in Fig. [2] show the transition energies for pure 28 Si, 
i.e. the transition energies between the unperturbed eigenstates obtained from our k • p 
calculations. For the maxima of the distributions calculated for nat Si, we observe a shift in 
energy with respect to the unperturbed case, which is significant for the outer Am=l and 
the Am=2 transitions. The maxima of the |1)-B-|2) and |3)-B-|4) distributions are shifted 
to higher energies by 0.315 /xeV and 0.378 /xeV with respect to 28 Si, respectively. For the 
|l)o|3) and |2)<B-|4) Am=2 transitions, we find shifts of 0.057/ieV and 0.151/ieV, respec- 
tively. These shifts manifest themselves in (temperature dependent) differences between the 
effective g-values of 28 Si and nat Si, which will be discussed further below. 

In order to see to which extent the properties of boron EPR spectra that have not been 
understood so far can be explained by the calculated distributions of the transition energies, 
we directly compare our results with experimental data. The open circles in Fig. |3(a) show a 
typical EPR spectrum of the Am=l transitions measured on a nat Si sample at 5| | [001]. The 
spectrum shows a broad structured line at 5=567.2 mT, which is attributed to a superposi- 
tion of the resonances originating from the |1)-H-|2) and the |3)-B-|4) transitions. 1,0 Further, 
we observe a narrower asymmetric resonance at 5=621.1 mT, which can be assigned to the 
|2)-B-|3) transition.— The origin of the narrow substructure feature centered at 5=570.2 mT 
and discussed under (iii) in Sec. His the result of an interplay between spin excitations and 
spin relaxations in the four level spin system and is the subject of a different publication.— 
In Fig. 11(a), the open circles show the EPR spectrum measured on the same nat Si sample in 
a magnetic field range where the Am=2 resonances are expected. The broad peak observed 
at 5=297.5 mT can be assigned to a superposition of the |1)-B-|3) and |2)-B-|4) resonances.— 
Here, the narrow substructure line has the opposite sign as compared to the substructure of 
the broad Am=l peak. The additional signal in the range from 5=314 mT to 5=324 mT 
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FIG. 3: (a) EPR spectrum of the Am=l transitions measured on a nat Si:B sample (open circles) 
shown together with the theoretical lineshape derived from the distribution of transition energies 
shown in Fig. [2] using A.E 29 =0.74 meV and A£' 30 =1.46 meV (solid line). The dashed curve 
shows the simulated spectrum for Ai? 29 =0.68 meV and Ai? 30 =1.34 meV. (b) Corresponding EPR 
spectrum of isotopically purified 28 Si:B. In (b), the solid line is a numerical fit using the first 
derivative of two Lorentz lines. The orientation of the external magnetic field is -B||[001] in both 
cases. The measurement temperature was 3 K. The dotted lines indicate the resonance fields (zero 
crossings) of the inner and outer Am=l resonances for nat Si. 

is not related to boron but is due to a background signal of our nat Si sample which is of 
unknown origin. 

The curves (b) in Figs. [3] and @] show the corresponding spectra that were measured on 
the 28 Si sample. We observe the same features as in the spectra of the nat Si sample, however, 
the linewidths and the line shapes of all resonances are strongly changed. The solid lines in 
Figs. |3jb) and show computer simulations of the spectra using single Lorentzian lines 
for each resonance peak, neglecting the narrow substructures. The best fit which leads to 
a very good agreement with the experimental data, is obtained for peak-to-peak linewidths 
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FIG. 4: (a) EPR spectrum of the Am=2 transitions measured on a nat Si:B sample (open circles) 
shown together with the theoretical lineshape derived from the distribution of transition energies 
recalculated for B=0.3 T using A£ 29 =0.74 meV and A£' 30 =1.46 meV. (b) Corresponding EPR 
spectrum measured on isotopically purified 28 Si:B. In (b), the solid line is a numerical fit using the 
first derivative of one Lorentz line. The orientation of the external magnetic field is £?||[001] in 
both cases. The measurement temperature was 3 K. The dotted line indicates the zero crossing of 
the Am=2 resonance for nat Si. 

of 7.2 mT and 1.7 mT for the low-field and the high- field Am=l resonances, respectively. 
For the Am=2 resonance, a linewidth of 3.1 mT is obtained. In spite of the Lorentzian 
line shapes, we conclude from the temperature dependence of the linewidths (not shown) 
that the |1)<H-|2), |3)f-)-|4) resonance and the Am=2 resonance shown are inhomogeneously 
broadened, since the linewidths of these two resonances become almost temperature inde- 
pendent below 3 K. We attribute this broadening to a distribution of the transition energies 
between the different energy levels induced by random local strains due to C, O, and B 
point defects.— 1 ^ A quantitative estimation of the point defect concentration is given later 
in the text. For the narrow Am=l resonance at higher magnetic fields, belonging to the 
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|2)-B-|3) transition, the strain induced broadening is only a second order perturbation. For 
this resonance, the linewidth continues to decrease for temperatures down to 0.3 K and at 
the temperature investigated here (3 K), the linewidth of this resonance is therefore a good 
estimate for the natural linewidth. 

Therefore, we use the line shape of the |2)-B-|3) resonance in 28 Si for a convolution with 
the statistical distribution of spin packets shown in Fig. [2] in order to obtain simulated line 
shapes for nat Si that can be compared to the experimental spectra of the Am=l and Am=2 
resonances in Fig. |3]^a) and Fig. HJ^a), respectively. To make a direct comparison of the line 
shapes easier, the calculated spectra have been shifted rigidly on the magnetic field axis 
so that they overlap with the experimental data. This shift corrects for the inaccuracy of 
the absolute g-values obtained from our calculations. However, we only performed a rigid 
shift of the entire spectrum by a constant field AB and the splitting of the Am=l reso- 
nances has not been changed. As our model does not yield EPR transition probabilities, 
the amplitudes of the simulated peaks have manually been adjusted to the measured values. 
Neglecting the substructure of the broad transitions caused by dynamic effects,— the calcu- 
lated signals (black solid curves) show a striking agreement with the experimental spectra, 
which demonstrates that all inhomogeneous broadening effects in our nat Si sample can be 
explained by isotope-induced fluctuations of the valence band edge alone. The apparent 
noise in the simulated spectra is the result of the still somewhat limited number of random 
isotope distributions used for the simulation. As indicated by the dotted vertical lines, also 
the shifts in the effective g-values between 28 Si and nat Si that have been predicted as shifts in 
the transition energies in Fig. [2} are found experimentally. For the broad Am=l resonance 
in Fig. |3l we observe a shift of 4.2 mT which translates into a shift in transition energy of 
0.3 /j,eV at 5=0.6 T and is in good agreement with the theoretical values. For the narrow 
Am=l resonance, we do not see a change in transition energy within the resolution of our 
measurement, as expected from our calculations. The Am=2 resonance shows a smaller shift 
of 0.5 mT corresponding to a change in transition energy of 0.1 fieV at 5=0.6 T. We note 
that the resonance fields of the substructures, both for the Am=l and Am=2 transitions, 
are not shifted due to the isotopic perturbation, which means that the substructure lines are 
not located exactly at the maximum of the broader absorption lines but are slightly shifted 
to higher fields in nat Si. 

We note that the narrow substructures on the broad Am=l resonance showing a negative 
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sign and the Am=2 resonance showing a positive sign for the temperature investigated here, 
do not originate from a superposition of two absorption peaks which are shifted with respect 
to each other on the magnetic field axis. Such a splitting of resonance lines could in principle 
result from slightly different g-values of the two boron isotopes 10 B and 11 B, which have a 
slightly different binding energy (1.9xl0~ 2 meV higher for 10 B).— However, due to their 
natural abundance of 19.9 % and 80.1 %, respectively, such an effect has to lead to an 
asymmetric resonance line, which contradicts the highly symmetric line shape of the B- 
related resonances observed in 28 Si, and can therefore be excluded. Another effect that 
in principle has to be considered for the understanding of the line shape is the hyperfine 
interaction of the acceptor holes with the nuclear spins of the corresponding acceptor nuclei 
and 29 Si ligands. For 10 B and 11 B, the nuclear spin is 3 and 3/2, respectively. If the 
acceptor hyperfine splitting could be resolved in our spectra, we would therefore expect a 
multiplet of lines rather than a twofold splitting of the resonance line that might be invoked 
to account for the central substructures. We neither observe a resolved hyperfine signature 
in any of our spectra nor do we need a further source of inhomogeneous broadening for the 
description of our spectra that could be attributed to unresolved B hyperfine multiplets or a 
strong superhyperfine interaction with 29 Si ligands. We conclude that the hyperfine coupling 
strength of B acceptors to their nuclei is much smaller than for the case of shallow donors 
in Si. This is a result of the p-character of valence band Bloch states, which cancels the 
strong Fermi contact hyperfine term at the position of the dopant atom due to a vanishing 
probability amplitude of the acceptor hole at its nucleus. For heavy hole states in III-V 
semiconductor quantum dots, it has recently been pointed out that the anisotropic dipolar 
hyperfine coupling strength can amount to up to ~10 % of the Fermi contact term of 
electrons.— However, it is not necessary to take into account an additional dipolar hyperfine 
or superhyperfine term for the description of our data. 

In Fig. [5](a), the FWHM of the sum of the transition energy distributions of the outer 
Am=l transitions (|1)<B-|2) and |3)-(->|4)) (circles), and the distribution of the |2)<B-|3) tran- 
sition energies (diamonds) are shown as a function of the degree of isotopic purification. 
As an approximation, the ratio between the residual 29 Si and 30 Si concentrations was kept 
constant at the nat Si value. For both resonances, we observe a monotonous decrease of the 
isotope- induced broadening when the 28 Si concentration is increased. In Fig. [5](b), the line 
shape of the distributions of the outer Am=l transition energies are exemplarily shown for 
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FIG. 5: (a) FWHM of the sum of the transition energy distributions of the |l)o|2) and |3)o|4) 
transitions (circles), and the distribution of the |2}o|3) transition energies (diamonds) as a function 
of the 28 Si concentration. The solid lines are guides to the eye. The calculations where performed 
for an external magnetic field B oriented along the [001] axis. The ratio between the residual 29 Si 
and 30 Si concentrations was kept constant at the value of nat Si. (b) Distribution of the |l)o|2), 
|3}o|4) transition energies for a 28 Si concentration of 99.2 % (dashed line) and 99.9 % (solid line). 
The dotted lines are numerical fits using a Gaussian line shape. A numerical fit to the 99.2 % 28 Si 
distribution using a Lorentzian line shape leads to an indistinguishable agreement with the actual 
distribution (solid line). 

99.2 % (dashed line) and 99.9 % (solid line) 28 Si. Both distributions are shown together 
with numerical fits assuming Gaussian lineshapes (dotted lines). For the 99.9 % 28 Si dis- 
tribution, a Lorentzian fit leads to a curve which cannot be distinguished from the actual 
distribution (solid line). Comparing the calculated distributions with the fits, it can be 
seen that the line shape of the distributions changes from Gaussian-like to Lorentzian in the 
regime of high isotopic purity with >99 % 28 Si. This complicates the discrimination between 
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isotope- induced and strain- induced broadening in the high 28 Si concentration range, where 
both types of inhomogeneous broadening are Lorentzian. For 99.98 % 28 Si, the expected 
FWHM due to isotope effects for the |2)-H-|3) transition and for the sum of the |1)-B-|2) and 
|3)o|4) distributions has decreased to 0.002 /ieV and 0.02 fieV, respectively. For the spec- 
trum shown in Fig. E(b), this means that the isotope- related inhomogeneous contributions 
to the observed linewidths are only 0.24 mT for the low-field resonance and 0.024 mT for 
the high-field line. These values fall more than one order of magnitude below the experi- 
mentally observed linewidths. The linewidth of the broad Am=l resonance in the spectra of 
our 28 Si sample should therefore only be determined by spin relaxation and the concentra- 
tion of point defects inducing random strain fields. Under this assumption, we can directly 



use Eq. 6 in Ref. [Ill to estimate the total concentration of C and O point defects from the 
peak-to-peak linewidth at low temperatures. We find a value of 3(±2)xl0 16 cm" 3 , which 
is consistent with the ^lxlO 16 cm -3 concentration of C impurities measured by infrared 
absorption spectroscopy on the ~607 cm -1 C local vibrational mode.— A further reduction 
of the total point defect concentration to ^lxlO 15 cm" 3 would be necessary to observe an 
isotope induced broadening in the 28 Si sample investigated. 

To obtain a more general picture, we have further investigated the anisotropy of the 
isotope-induced broadening effects. We have calculated the broadening of the Am=l EPR 
lines for different directions of the magnetic field ranging from S||[001] (6=0°) to S||[110] 
(0=90°). In Fig. [6(a), the peak-to-peak linewidth of the broad Am=l resonance (|1)<B-|2) 
and |3)-B-|4)) calculated for nat Si using AE 29 =0.74 meV and AE 30 =1.46 meV is shown as 
a function of the angle 6 (solid cirlces). We observe a strongly anisotropic behavior with 
a minimum of the linewidth at _B||[001] and a maximum at S||[lll]. A comparison with 



experimental data from Ref. 



11 



(open squares) and those obtained in this work (open tri- 
angles) shows a very good agreement with our theoretical calculations. The experimental 
uncertainty should be the smallest for small values of 9 where narrow lines with a large 
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show a corn- 



amplitude are observed. Regarding this, the experimental data from Ref. 
paratively large deviation from our theoretical result for 0=0° and #=20°. These larger 
linewidths can however be understood from a higher concentration of C and O point defects 
in the sample measured by Neubrand, which imposes a lower limit to the experimentally 
observed linewidth via an additional inhomogeneous broadening induced by random local 
strain. 
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FIG. 6: Linewidth anisotropy of the |1}-H-|2), |3}-o-|4) EPR transitions (a), and the |2)-f->-|3) transi- 
tion, shown in (b). The calculated angular dependence (full circles) for nat Si using A.E 29 =0.74 meV 
and A£' 30 =1.46 meV is compared with experimentally determined linewidths for nat Si from this 
work (triangles) and the anisotropy obtained from Fig. 4 in Ref. llOl (squares). The lines are guides 
to the eye. 

In Fig. |6](b), the theoretical and experimental data obtained for the |2)-H-|3) resonance 
are shown. Compared to Fig. E](a), the anisotropy is inverted. We find a maximum of 
the linewidth for 5||[001] and a minimum for £?||[111]. Again, the experimentally observed 
linewidth anisotropy for nat Si is very nicely reproduced by our theoretical model. The origin 
of the relatively large scatter of the experimental data points for 6<30° is the asymmetry 
of the |2)-B-|3) resonance shown in Fig. [3]^a). For S||[001], the positive amplitude a (with 
respect to the zero crossing) on the low-field side of the derivative line is smaller than the 
negative amplitude b. This asymmetry of the line shape, which results from the mixing of 
the Zeeman states due to the isotopic perturbation, also leads to a relatively weak and broad 
negative peak on the high field side, which, for a given noise level, increases the experimental 
error for the determination of the overall linewidth. In Fig. [6](b), we find the best agreement 
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FIG. 7: Angular dependence of the asymmetry ratio a/b of the 1 2) -<— >■ 1 3) transition as defined in 
Fig. [3l The triangles show the experimentally obtained values for the nat Si sample studied here, 
full circles represent the theoretical results. The solid line is a guide to the eye. 

between theory and experiment around 1 1 [111], where we observe the smallest linewidths. 

As shown in Fig. [TJ also the asymmetry a/b of the |2)-B-|3) resonance line shape has a 
minimum at this orientation. Triangles show the experimental data, full circles represent the 
asymmetry that is deduced from the simulated line shapes. With the exception of two data 
points around #=70°, the angular dependence of the line shape anisotropy can also quan- 
titatively be understood from our theoretical model. As mentioned above, the asymmetric 
broadening of the 1 2) -<— >• 1 3) resonance originates from a mixing of the unperturbed Zeeman 
states j^j) via the off-diagonal matrix elements of the isotopic perturbation Hamiltonian 
H^ 4 . The difference between the experimental and the theoretical values of the asymmetry 
factor in the region of 0=70° could be a result of a slight anisotropy of the linewidth of the 
single spin packets contributing to the inhomogeneously broadened line, which has not been 
taken into account for the simulation. The scatter around a/b=l in Fig. [7] lies within the 
experimental error bars resulting from the different spectral overlap of the resonance lines 
for different orientations of the magnetic field. 

The calculated distributions of energy splittings in Fig. [2] and the experimental EPR 
spectra in Figs. [3] and Fig. H] revealed that the maxima of the distributions of transition 
energies in nat Si are shifted with respect to the corresponding transition energies observed 
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FIG. 8: (a) Theoretical anisotropy of the effective g-values of the outer Am=l transitions as 
obtained from our k • p model for 28 Si (solid symbols) and including the isotopic perturbation in 
nat Si (open symbols). Diamonds and triangles indicate the g-values obtained for the |l)o|2) and 
|3)-f->|4) transitions, respectively. The solid and dashed lines are numerical fits to the data using 



the theoretical anisotropy given by Eq. 5 in Ref. 



1CJ as a model. In (b) and (c), a direct comparison 



of the experimental and the theoretical g-value anisotropics for both Am=l resonances is shown. 
Data points obtained for 28 Si are indicated by open circles, nat Si data points are shown as open 
triangles. The theoretical data for the outer Am=l resonance is a weighted average of the |l)o|2) 
and |3)-f-)-|4} transitions shown in (a) for T=3 K. The crosses show the g-values of the substructure 
line. Details concerning the calculations of the data shown in (c) are described in the text. 

for 28 Si. In Fig. [8]^a), we show how this effect depends on the orientation of the magnetic 
field for the outer Am=l transitions. The effective g-values calculated from the theoretical 
|1)<H-|2) and |3)-H-|4) transition energies in pure 28 Si are shown as full diamonds and trian- 
gles, respectively. The open symbols indicate the corresponding effective g-values for nat Si 
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calculated from the maxima of the transition energy distributions. According to Ref. 
assume the g-values of B in Si to be negative. As guides to the eye, we also show numerical 
fits to the data based on the theoretical anisotropy given by Eq. 5 in Ref. |l0| as solid and 
dashed lines. For 28 Si, we observe a splitting in energy of the two outer Am—1 transitions 
for 9 > 0°. This splitting has a maximum for 5||[111] and should result in a splitting of the 
broad Am=l EPR resonance. However, we cannot resolve this splitting in our experiments 
due to the comparatively strong defect-induced inhomogeneous broadening. For nat Si, the 
g-values are shifted to more negative values with respect to 28 Si for both transitions. The 
|3)-f->-|4) resonance is shifted more strongly than the |l)-f^|2) resonance. 

In our EPR experiments we observe only one resonance with a g-value that corresponds 
to an average of the g-values of the two outer Am=l transitions that is weighted according 
to the equilibrium populations of the different spin states. In Fig. Mb), the experimentally 
obtained g-value anisotropies of the outer- as well as the inner Am=l resonances for 28 Si 
(open circles) and nat Si (open triangles) are plotted. In addition, also the effective g-values 
of the substructure superimposing the outer Am=l resonance are shown for nat Si (crosses). 
For direct comparison, the corresponding theoretical results are shown in Fig. EJc). The 
effective g-values of the outer Am—1 resonances were calculated by taking the maxima 
of the weighted sums of the |1)-H-|2) and |3)-H-|4) transition energy distributions assuming 
Boltzmann distributed populations of the four spin states (T=3 K). The effective g-values 
at which the substructure line on the broad Am—1 resonance is expected were calculated 
under the assumption that the substructure lines result from a subensemble of acceptors for 
which the transition energies of the two outer Am—1 transitions are equal, i.e. the energy 
level scheme of the four level spin system is symmetric— The solid lines in Figs. Mb) an d 
(c) are again numerical fits using the model described in Ref. 1(1 We note that the g- 
value scale in Fig. [S^c) is shifted by 0.04 with respect to the one in Fig. Mb) to enable a 
better comparison between the experimental and theoretical data. Disregarding this rigid 
shift, which is caused by the error of the absolute g-values obtained from our k ■ p model, our 
calculations can be used to describe and understand the different experimental observations. 
For the inner Am—1 transition |2)<B-|3), we could already see in Figs.[2]and[3]that the isotopic 
perturbation does not lead to a measurable shift in the g-value for the 0° orientation. As 
can be seen in Fig. M this finding holds for all orientations of the magnetic field. Comparing 
the anisotropies of the outer Am=l resonance in Figs. Mb) and (c) for 28 Si and nat Si, we 
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observe a similar shift of the g-values. Due to the large linewidth of the latter resonance of 
up to 100 mT in nat Si, it was not possible to extract a reliable g- value for all intermediate 
orientations of the magnetic field from our experimental spectra. 

The g-values of the substructure line in nat Si, which could be determined for all orien- 
tations of the magnetic field and with a higher precision due to its comparatively small 
linewidth, show only a small deviation from the effective g-values of the outer Am=l reso- 
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this 



nance in 28 Si irrespective of 9. In clear contrast to what has been reported in Ref. 
means that the resonance field of the substructure line does not coincide with the center 
of the broad Am=l resonance in nat Si. The agreement between our theoretical model and 
the experimental data in this point strongly indicates that the substructure lines indeed 
originate from a subensemble of acceptors where the energy splittings between the states |1) 
and |2) and the states |3) and |4) are equal. 

For the Am=l substructure line in nat Si, the best fits to our experimental data is obtained 
for gi=-1.0776 and q2=-0-0307. For the inner Am=l resonance, we obtain pi =-1.0728 and 
(72=-0.0315. In Ref. lid, for comparison, pi =-1.0740 and (?2=-0.0307 was determined for the 
Am=l substructure and gi=-1.0676 and g2=-0.0317 was measured for the inner Am=l 
resonance. We find a good agreement of our data with the original reports in particular for 
g> 2 , which determines the angular dependence of the effective g-values. For gi, which causes 
an angle-independent offset of the g-value anisotropy, we obtain slightly lower values. This 
deviation might result from differences in the calibration of the magnetic field. 

Figures |3] to [7] show that the line shape and width of the Am=l and Am=2 resonances 
in nat Si with a low concentration of point defects can quantitatively be understood by our 
model accounting for the fluctuations of the valence band position caused by the disorder in 
the placement of the different Si isotopes in the crystal lattice in the vicinity of the different 
B acceptor nuclei. In contrast, the original model used to calculate the acceptor ground state 
splitting for B=0 could not explain the experimentally observed value quantitatively— In 
Fig. El the calculated statistical distribution of the B acceptor ground state splitting taken 



from Ref. 



14 



(dashed line) is shown together with experimental data that were obtained from 
phonon absorption spectroscopy.— 1 ^. The shapes of the distributions are in good agreement, 
however, the maximum of the splitting energy shows a discrepancy of 30 %. To investigate 
the origin of this discrepancy, we have repeated the calculations for the distribution of the 
residual ground state splitting in nat Si applying our theoretical method to B=0. The result 
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FIG. 9: Energy distribution of the residual B acceptor ground state splitting in the absence of 
an external magnetic field. The open circles show the combination of two sets of experimental 
data measured with Si samples with a boron concentration of 5.4xl0 15 cm" 3 and 2xl0 14 cm"^ 
The data is taken from Fig. 3 in Ref. 



24l . The dashed line represents the distribution calculated 



in Ref. |14| . The solid line shows the results obtained with our theoretical method using the same 
values for the valence band offsets between the different Si isotopes as in Ref. and the dotted 
line is obtained when AE 29 and AS 30 are decreased by 8 % 

is depicted as the solid line in Fig. [9] and shows a much better quantitative agreement with 
the experimental data. Besides the potential Q at the central grid node, which was not 
adjusted for the calculations for B=0, the most relevant parameters in our model are the 
valence band offsets between the different Si isotopes (Vi SO ). Since the values used for these 
parameters in this work so far are identical to the ones used in Ref. [14], the discrepancy 
between the acceptor ground state splitting at B=0 calculated in the original work and in 
the present_work can not be attributed to an uncertainty in the valence band offsets deduced 
i It should rather result from a less realistic modeling of the boron acceptor 



from Ref. 
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wave function. In order to see how sensitive the residual ground state splitting depends on 
the assumed valence band offsets, we have recalculated the distribution of splitting energies 
for different AE 29 . The second parameter AE 30 is directly linked to the value of AE 29 by 
the M _1//2 dependence of the renormalization energy on the isotope mass M. Using our 
wave functions, the best overall agreement with the experimental data shown in Fig. |9] is 
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obtained when AE 29 and AE 30 are decreased by 8 % from their original value (dotted line). 
We have also recalculated the EPR spectrum of the Am—1 resonances in nat Si with the 
reduced AE values and the result of this simulation is shown as the dashed line in Fig. [3](a). 
Using the lower band offsets, we obtain a slightly worse agreement with the experimental 
spectrum on the low-field edge of the broad |l)-f-)-|2), |3)-B-|4) resonance, while an equally 
good or even slightly better agreement is found for the central region and the high field 
edge of this resonance as well as for the |2)<B-|3) line. In view of the overall agreement with 
the experimental data and neglecting possible inaccuracies of the wave functions obtained 
from our calculations, we suggest a slight correction of the valence band offsets between the 
different Si isotopes to A£' 29 =0.68 meV and Ai? 29 =1.34 meV. We note that a deviation 
from the original values by more than 5 % to higher energies and more than 10 % to lower 
energies would lead to a significant discrepancy between theory and experiment. 



V. CONCLUSIONS 



We have set up a theoretical model that allows for a quantitative description of the iso- 
topic effect on the Zeeman energy level scheme of paramagnetic B acceptors in Si, accounting 
for the isotopic perturbation via local fluctuations of the valence band edge at 29 Si and 30 Si 
lattice sites. Using this model, we can reproduce the experimentally observed line shapes of 
the B-related Am—1 and Am=2 resonances in nat Si with excellent agreement between theory 
and experiment. Our results demonstrate that the inhomogeneous broadening of B-related 
EPR lines measured in ultra-pure natural silicon can quantitatively be understood from the 
isotopic perturbation, both in its magnitude as well as in the resulting line shapes, without 
the necessity to invoke hyperfine interaction with the B nucleus or with 29 Si. Concretely, the 
open questions raised in the introduction can be answered as follows: (i) As can directly be 
deduced from a comparison of the data shown in Fig. [6] and Fig. 2 in Ref. |UJ, the Gaussian 
contribution to the line shape of the broad Am=l resonance is merely a phenomenological 
way to account for the broadening induced by the isotopic randomness. However, a Gaussian 
line is no realistic representation of the isotope-induced broadening, which is determined by 
the distribution of transition energies shown in Fig. [2l Therefore, an additional Lorentz line 
was needed in Ref. lllj to obtain a reasonable fit of the broad Am=l resonance measured 
with the samples of highest crystalline purity. Consequently, the 10 mT threshold of the 
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Lorentzian contribution to the line shape (ii) was an artifact that resulted from an inappro- 
priate model of the line shape in the limit of low point defect concentrations. We have found 
that the perturbation due to the random distribution of the different Si isotopes also leads 
to a shift in the effective g- values, mainly of the broad Am=l and the Am=2 resonances. 
We note that this isotope induced shift could in principle be used to control the g- values of a 
B acceptor by tailoring its isotopic surrounding. This could also be achieved dynamically by 
placing an acceptor close to an isotopic heterojunction and manipulating its wave function 
with the help of electric fields. This idea might particularly be of interest for applications 
e.g. in quantum computation technology, as it could be realized in a completely nuclear spin 
free environment only using 28 Si and 30 Si isotopes. To give an example, we obtain an average 
g- value of the broad Am=l resonance of -1.267 for an isotopically engineered material com- 
posed of 50 % 28 Si and 50 % 30 Si (S||[100], T=3 K), which corresponds to a 4.2 % change 
of the g-value compared to g=-1.216 calculated for pure 28 Si. Concerning issue (iii), our 
results suggest that the substructure lines of the Am=l and Am=2 resonances originate 
from a subensemble of acceptors, where the Zeeman level scheme is perturbed such that the 
|1)-H-|2), |3)<H-|4) transition ener gies are equal. We discuss these substructure lines and their 



dynamics in more detail in Ref. 
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Finally, a comparison of our calculation with previous 
works investigating the B acceptor ground state in the absence of an external magnetic field, 
provides an independent verification of the energy offsets between the valence bands of the 
different isotopes of silicon. 
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